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ABSTRACT 

HD molecules can be an important cooling agent of the primordial gas behind the 
shock waves originated through mergings of the dark matter haloes at epochs when 
first luminous objects were to form. We study the necessary conditions for the HD 
cooling to switch on in the low temperature range T < 200 K. We show that these 
conditions are fulfiled in merging haloes with the total (dark matter and baryon) mass 
in excess of M cr ~ 10 7 [(1 + z)/20]~ 2 M & . Haloes with masses M > M cr may be the 
sites of low-mass star formation. 

Key words: early Universe - cosmology:theory - shock waves. 



1 INTRODUCTION 

Cooling of primordial gas is determined by molecular hydro- 
gen, which forms in an expanding universe after recombina- 
tion (Lepp & Shull 1983, Puy et al 1993, Palla, Galli & Silk 
1995, Galli & Palla 1998, Stancil, Lepp & Dalgarno 1998). 
Inside the first virialized dark matter haloes the fraction of 
H2 can reach 10 -3 , which is able to cool gas to 200 K. At 
lower temperatures energy losses in roto- vibrational H2 lines 
become insufficient to cool gas further. On the other hand, 
owing to a low rotational energy and a large dipole moment 
HD molecules might be an efficient coolant at T < 200 K. 
However, the abundance of HD is fairly sensitive to ther- 
mal history of the gas, so that the question of whether HD 
dominates radiative cooling at low temperatures or its con- 
tribution is negligible depends on physical conditions along 
the evolutionary path. This circumstance is reflected in con- 
radictory conclusions about the role of HD in thermody- 
namics of primordial gas: for instance, Lepp & Shull (1983), 
Bromm, Coppi & Larson (2002) point out that HD cooling 
cannot dominate in primordial gas when the first objects 
form, while Bougleux & Galli (1997), Puy & Signore (1997, 
1998), and more recently Uehara & Inutsuka (2000), Flower 
(2002), Nakamura & Umemura (2002), Flower & Pineau 
des Forets (2003) and Machida et al (2005) find that HD 
molecules form in significant amount and can play a domi- 
nant role in thermal evolution of pregalactic gas. 

Transformation of H2 molecules into their isotop ana- 
log HD is energetically favoured because of higher binding 
energy of HD molecules. In equilibrium (Solomon & Woolf 
1973, Varshalovich & Khersonskii 1976) 
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where factor 2 stems from the difference between the chem- 
ical constants of HD and H2 molecules: x H d ~~ Xh 2 ~ m ^ 
(Landau & Lifshitz 1969). It is seen, that at T 150 K ra- 
diative HD cooling can enhance chemical fractionation of 
deuterated molecules progressively: a small decrease in tem- 
perature results in an increase of n(HD)/n(H2), what in turn 
increases radiative cooling in HD lines (Shchekinov, 1986). 
In realistic conditions thermal equilibrium can be reached 
only in a characteristic time depending on a set of chemical 
reactions which control the transformation H2 ^=>HD. In 
primordial gas H2 converges into HD most efficiently in the 
reaction H 2 +D+ -> HD+ H + (Palla et al. 1995); other chan- 
nels are unimportant on the Hubble time. The background 
H2 abundance and the frozen-out post-recombination frac- 
tional ionization x ~ 10~ 4 are sufficient to support deuter- 
ated H2 only at the level n(HD) 10 _3 n(D). Thus, a neces- 
sary condition for the H2 — > HD convergence to be efficient 
is that both H2 and electrons are abundant. In a cosmologi- 
cal pregalactic substrate such conditions are most naturally 
fulfiled behind shock waves associated with formation and 
virialization of dark matter haloes (Tegmark et al. 1997, 
Barkana & Loeb 2001). When formed in sufficient amount 
HD molecules can cool gas down to the temperature of the 
CMB radiation - the lowest possible temperature, and thus 
provide conditions for formation of low-mass stars. It is clear 
that this can occur with shock waves whose parameters (for 
instance, postshock temperature, total mass involved, and 
so on) have reached some critical values. 

In the hierarchical scenario dark matter haloes of larger 
mass form through mergings of smaller objects - minihaloes 
(see for review Barkana & Loeb 2001, Ciardi & Ferrara 
2004). In these conditions the baryonic component of dark 
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haloes is compressed by shock waves. Apparently, all subse- 
quent evolution of baryons in dark haloes and their conden- 
sation into stars is determined by these shock waves. Futher- 
more, one can consider that the virial state of dark matter 
haloes and their baryons is reached through formation of 
shocks. It is clear, that the physical state of the baryons 
processed by shocks depends on many factors, such as ge- 
ometry of a shock front, duration of a compressing flow etc. 
However, basics characteristic features can be understood in 
the framework of a head-on collision of two clouds of equal 
masses (Suchkov, Shchekinov & Edelman 1983, Shapiro & 
Kang 1987). In the context of the hierarchical scenario a 
possible role of shock waves in enhancing H2 cooling and 
stimulation of primordial star formation was studied by Ya- 
mada & Nishi (1998), and more recently by Cen (2005). 
Radiative cooling in H2 lines was indeed found extremly ef- 
ficient in decreasing the postshock temperature to ~ 100 
K. In spite of the fact that HD molecules are expected to 
form at this temperature, HD chemistry was not taken into 
consideration in this analysis. The main goal of our study 
are the conditions when HD molecules can be the dominant 
cooling agent in a pregalactic gas. 

In Section 2 we describe a simple thermo-chemical 
model of a shocked gas; in Section 3 the conditions for HD 
molecules to be the dominant cooling agent are analyzed; 
in Section 4 we discuss qualitatively a possibility of a cold 
compressed layer to fragment; the discussion and summary 
of the results are given in Section 5. 

Throughout the paper we assume a ACDM cos- 
mology with the parameters (f2o, £Ia, f2m, fif>, h) = 
(1.0, 0.71, 0.29, 0.047, 0.72) as inferred from the Wilkin- 
son Microwave anisotropy Probe (WMAP), and deuterium 
abundance 2.62 x 10~ 5 (Spergel et al. 2003). 



2 THERMO-CHEMISTRY OF POSTSHOCK 
GAS 

For simplicity, we assume a head-on collision of two identi- 
cal minihaloes, which are virialized at redshift z and have 
density of matter 187r 2 f2 m/ o c (l + z) 3 , p c being the critical 
density; the collision is assumed to take place at the same 
redshift z. The minihaloes move with the relative velocity 
v c = v3<7, where a is a ID velocity dispersion of the larger 
halo 

° 2 = Sr = GM 2/3 (3TT :i n mPc ) 1/3 (i + z). (2) 

In the center of mass a discontinuity forms at the symmetry 
plane, and two shocks propagate outwards. The time for the 
shock to pass through the entire minihalo is t c = 3R/2v c . As 
it was pointed out by Gilden (1984) for supersonic collisions 
the rarefaction time in transverse direction is greater than 
the collision time, which means that a ID scheme describes 
the overal picture qualitatively correct. We assume the post- 
shock flow to be isobaric (Suchkov et al. 1983, Shapiro & 
Kang 1987, Anninos & Norman 1996). The initial post-shock 
temperature is 



where the post-shock gas is assumed to be isothermal, i.e. 
temperatures of the electrons, ions and neutrals are kept 
always equal T e — Ti = T„. 

Chemical kinetics and thermal evolution of the gas are 
described by the following set of equations 

Xi = Fi(x,T,n)-Di(x,T,n), (4) 
f= ^^(r i (x,T,n)-A l (x,T,n)) + ^n, (5) 

i 

where k is the Boltzmann constant, Xi is the fraction of i-th 
species, Fi(x,T, n), Di(x,T, n) are respectively the forma- 
tion and destruction rates, Fi(x, T, n), Aj(x, T, n)) are the 
heating and cooling rates; x includes 14 species: H, H + , H _ , 
He, He+, He++, H 2 , H+, D, D+, D", HD, HD+, e. In the 
energy equation all relevant cooling processes (such as 
collisional excitation and ionization, Compton cooling, re- 
combination and bremsstrahlung radiation, cooling by HI, 
H2 and HD line emission) are accounted for. For the cooling 
functions of H2 and HD we adopted the expressions given 
by Hollenbach & McKee (1979), and Flower et al. (2000), 
respectively; note, that the Hollenbach & McKee H2 cooling 
function practically coincides with the one calculated by le 
Bourlot, Pineau des Forets & Flower (1999). The most recent 
calculations of the HD cooling function by Lipovka, Nunez- 
Lopez, & Avila-Reese (2005) show that in the temperature 
range of interest (T < 10 3 K) it coincides with that found by 
Flower et al. (2000). At higher temperatures, where Lipovka 
et al. (2005) predict an order of magnitude enhanced cooling 
from HD, the abundance of HD molecules is too low to con- 
tribute significantly. The other cooling functions connected 
with excitation and ionization of H and He are taken from 
Cen (1992). In general, radiation from a hot shocked gas can 
heat the gas, however, for the shock waves with the veloc- 
ity < 60 km s" 1 (Shull & McKee 1979, Shull & Silk 1979) 
its contribution to ionization and heating is unimportant, 
therefore in our calculations we assume T;(x, T, n) — 0. The 
chemical reaction rates are taken from Galli & Palla (1998). 
The thermo-chemical equations 0, @ are solved in one 
collisional time t c ; t c is much shorter than the age of the 
universe at any redshift z. 

Fig.l shows an example of physical characteristics of 
the postshock gas versus temperature in the collision of two 
minihaloes with the total masses M — 10 7 M© each at 
z — 20. The fractional ionization x and H2 abundance be- 
hind the shock are taken equal to their background values 
x = 10~ 4 and /(H 2 ) = 10~ 5 , /(HD) = 10~ 9 . Immediately 
behind the shock concentration of H2 increases very rapidly 
and reaches its maximum value ~ 2 x 10 -3 promoting cool- 
ing of gas to T ~ 200 K. If H2 cooling decreases the gas 
temperature below T = 150 — 160 K the abundance of HD 
molecules grows fast and HD cooling becomes dominant. It 
sets a runaway cooling on. At T < 100 K almost all the 
deuterium is converted into molecular form, and in one col- 
lisional time gas cools down to T ~ 60 K. 

Haiman, Rees, & Loeb (1997) showed that the very 
first stars emit a large amount of photons in Lyman- Werner 
band (LW), which produces a hostile UV background for 
HD molecules. The characteristic time of photodissociation 
of HD can be estimated as t^iss = 10 4 /Jlw yr, where Jlw 
is the LW flux in units of 10 -21 ergs s _1 cm -2 Hz -1 sr" 1 
(Le Petit, Roueff, & Le Bourlot, 2002, Johnson & Bromm, 
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Figure 1. Thermo-chemical evolutionary state vs temperature of a postshock gas parcell in the collision of minihalocs with the mass 
M = 10 7 Mq each at z = 20: in the upper panel abundances of the species: electrons (solid), H2 (dashed) and HD (dotted) lines - are 
shown; the lower panel shows the relative contribution of H2 (solid) and HD (dashed) molecules to the total cooling rate. 



2005). One can assume that the LW flux is by an order 
of magnitude equal to the ionizig flux from the first stars 
Jlw ^ J21 ~ exp[— (z — 5)] in the range 5 < z < 20 (Ciardi 
et al., 2000). On the other hand, the characteristic forma- 
tion time of HD is f form ~ 3 x 10 6 yr at z = 10 - 20. Thus, 
practically in the whole range of redshifts considered here 
photodissociation from the very first stars can be neglected: 
tdiss becomes comparable to £f rm only at z £ 10. At lower z 
self-shielding seems to come into play. Indeed, for the masses 
of merging haloes M > 1O 6 M0 the characteristic column 
density of HD molecules in the shocked layer at final stages 
can be as high as ~ 3 x 10 15 cm~ 2 , which is factor of 3 
greater than the column density needed to provide optical 



depth in LW band of the order of one < 
et al., 2002). 



10 15 cm" 2 (Le Petit 



3 HD BEHIND SHOCK WAVES 



Haloes of small masses (M < 5 x 10 6 M©) merge with 
slow relative velocities, and the post-shock temperature is 
so low that fractional ionization remains equal to its post- 
recombination value x ~ 10~ 4 . As a result, {x,T) dia- 
grams connecting x(t) and T(t) along evolutionary paths 
very weakly depend on halo masses, as seen from Fig. 2. 
Mergings of more massive haloes produce shocks with tem- 
perature T ^ 10 4 K, so that fractional ionization can be as 
high as x ~ 10~ 2 . Immediately behind the shock x starts 
growing from its pre-shock value to the corresponding equi- 
librium. However, as the characteristic collisional ionization 
time is longer than the cooling time, gradual increase of frac- 
tional ionization is accompanied by cooling, which in turn 
weakens the rate of collisional ionization. As a result, frac- 
tional ionization remains always lower that the equilibrium 
value determined by collisional ionization and radiative re- 
combination. At temperature T ~ 10 4 K x reaches the max- 
imum which depends on the postshock temperature, and 
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then goes down as shown in Fig. 2. This behavior restricts 
the efficiency of molecule formation: as seen from Fig. 1 H2 
saturates at T 10 3 K when fractional ionization begins 
to decrease sharply. HD instead monotonously grows until 
practically all deuterium is converted into molecular form. 
The abundance of H2 molecules at the end of baryon com- 
pression increases with the mass of the merging haloes, and 
the minimum temperature reached at this stage becomes 
lower. When it reaches T ^ 150 K a rapid fractionation of 
HD molecules sets on. For larger halo masses HD cooling 
dominates and cools baryons down to the temperature of 
the CMB radiation. Fig. 3 connects variations of tempera- 
ture and molecular concentrations (H2 and HD) along the 
evolutionary paths for different masses of the haloes merg- 
ing at z = 20. Their contribution to the total cooling is 
shown in Fig. 4 where open circles show the ending stage 
after one collision time t c : while immediately after cross- 
ing the shock front the contribution from HD (left curve) 
is negligible, it increases in the course of gas cooling, and 
the transition to the regime with a predominance of HD 
cooling is clearly seen at M ~ 10 7 M Q . The critical mass 
corresponding to the onset of a predominantly HD cooling 
is M™ ~8x 10 6 [(1 + 2)/2O]~ 2 M as shown in Fig. 5. This 
value is factor of 5 higher than the minimum mass of the 
first objects found by Tegmark et al. (1997). One may con- 
clude therefore, that at the lower end of masses the primeval 
gravitationally bound baryonic condensations stop cooling 
at relatively high temperature T ~ 200 K due to H2 radi- 
ation, and HD molecules do not form in sufficient amount 
to cool the baryons further. However, soon later, at z ~ 18 
more massive haloes can form and merge with a clear pre- 
dominance of HD cooling, and in principle with lower masses 
of possible stars. Note, that the critical mass M™ is fairly 
unsensitive to the choice of cooling function. In particular, 
for the cooling function of Galli & Palla (1998), who in- 
cluded the collisional rates from Martin, Schwarz & Mandy 
(1996) and Forrey et al. (1997), which in the temperature 
range T 10 4 K is approximately factor of 3 higher than 
H2-cooling function of Hollenbach & McKee in our calcula- 
tions, M™ increases only by 10-15% in comparison with 
that shown in Fig. 5. 



4 FRAGMENTATION OF THE COLD GAS 
LAYER 

When cooled the shocked gas layer can fragment due to grav- 
itational instability (Stone 1970, Elmegreen & Elmegreen 
1978, Gilden 1984). For a layer formed in a collision of two 
clouds the necessary conditions for instability imply that i) 
the characteristic growth time is shorter than the collision 
time, and ii) the critical wavelength is shorter than the ini- 
tial size of the clouds (Gilden 1984). The dotted line in Fig. 5 
depicts the critical mass of a merging halo gravitation- 
ally unstable in Gilden sense: merging haloes with individual 
masses above this line are unstable againts fragmentation; 
M° is approximated as M° ~ 10 7 [(1 + z)/20] _2 M Q . It is 
readily seen that M„ is practically coincident with the mass 
corresponding to the transition to a predominance of HD 
cooling M™ . The explanation of this coincidence can be 
found in the fact that the shock compressed layers become 
unstable in Gilden sense only when HD molecules cool the 



compressed gas below T ~ 150 K (Vasiliev & Shchekinov 
2005). It is obvious therefore that all masses M > at 
the latest stages are dominated by HD cooling, and as a 
consequence the fragments compressed by mutual action of 
ram pressure and gravitation are able to cool further down 
to T < 100 K. In the layer formed in a collision of two haloes 
with masses M = 10 7 M Q the corresponding Jeans mass is 
Mj = 10 3 M Q , which is less than 1% of the baryonic mass 
involved into collision. Fig. 6 shows the dependence of Jeans 
mass in unstable layers on the mass of merging haloes and 
the redshift z when the haloes merge: thick dotted white 
line delimits the haloes whose compressed layers are gravi- 
tationally unstable - all masses above this line are unstable. 
The thin solid lines are the iso- Jeans masses in compressed 
layers. When temperature in the compressed layer goes to 
the CMB value the Jeans mass becomes independent on red- 
shift. Clarke & Bromm (2003) have estimated a lower limit 
of the masses of fragments formed in similar conditions, with 
baryons confined by ram pressure from shocks and gravita- 
tion. However, for the primordial gas they assumed the lower 
temperature limit T = 200 K determined by H2 cooling. As a 
result, the masses of the fragments are systematically higher 
compared to our estimates. In order to obtain smaller values 
of the masses they had to assume that the compressed gas 
has cooled down to the CMB temperature due to cooling 
provided by CO molecules only after the gas is pre-enriched 
with metals. 

The lower panel of Fig. 6 shows the gas temperature 
at the latest stages of gas compression, i.e. at t = t c , for 
mergings of haloes with masses M at redshift z. It is clearly 
seen that the compressed gas can cool down to the lowest 
temperature Tcmb(z) in haloes with masses M > M™ as 
their cooling is dominated by HD molecules. Note, that this 
result does not contradict the conclusions drawn from nu- 
merical simulations by Abel et al. (2000, 2002) and Bromm 
et al. (2002), because these simulations have been done for 
halo masses close to the minimum mass expected for first 
objects. 

Fragmentation of a compressed gas layer can be ini- 
tiated by thermal instability at the stages when H2 and 
HD begin to dominate radiative cooling. In Fig. 7 (upper 
panel) temperature dependence of the effective ccoling func- 
tion A(T) and its logarithmic derivative $ = (d In A/d In T) n 
are shown in the temperature interval T < 10 4 K, where 
HI radiative energy losses become unimportant and molec- 
ular cooling sets on. The curves corresponding to different 
masses are plotted in temperature interval from the initial 
value immediately behind the shock to the lowest value at 
t = t c - In a non-steady cooling medium without external 
heating the condition for the isobaric mode to be unstable 
reads as <3> < 2 (e.g., Shchekinov 1978). It is evident that 
this condition is fulfiled at T £ 9000 K due to cooling of H2 
molecules, and in the lower temperature end at T <, 100 K 
where HD cooling dominates. The characteristic length of 
thermally unstable perturbations is Ar, ~ c s tr, where 



is the characteristic radiative cooling time, c a is the sound 
speed. The corresponding masses Mr ~ pr(Ar/2) , esti- 
mated for the parameters of a shocked gas, are shown in the 
lower panel of Fig. 7; here pn is the baryon density at the 
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Figure 2. The evolutionary paths connecting fractional ionization and gas temperature behind shocks in collisions of minihaloes of 
masses 10 6 , 2 X 10 6 , 3 X 10 6 , 5 X 10 6 , 8 X 10 6 , 10 7 (from left to right) at z = 20. The arrows show the direction of evolution of a gas 
parcel starting from its state immediately behind the shock; paths for the first four masses practically coinside. 




Figure 3. The evolutionary paths connecting variations of temperature and concentrations of H2 and HD molecules for the same masses 
as in Fig. 2 (from left to right). 
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Figure 4. The relative contribution of H2 (right curve) and HD (left curve) into the total cooling for the same set of models shown in 
Fig. 2: increasing numbers correspond to the increasing masses in Fig. 2; open circles mark the state after one collision time t c . 




Figure 5. The mass-redshift diagram delimitating interval of halo masses, where HD molecules dominate in cooling at low temperatures 
(this is the range above the dash-dotted line), and halo masses with the postshock baryon layer being gravitationally unstable and 
subjected to further fragmentation (these masses above the dotted line). The thick solid line shows the 3cr fluctuations in the ACDM 
cosmology; the dashed line shows the minimum mass limit found by Tegmark et al (1997). 



stage when thermal instability begins. For the haloes of large 
masses the mass of fragments which can form through ther- 
mal instability is lower than the total mass of baryons and 
Jeans mass practically in the whole temperature range of 
compressed gas layers. This means that in these conditions 
thermal instability stimulates subsequent fragmentation of 



the compressed gas. The minimum mass of fragments ex- 
pected from thermal instability Mr ~ 1O 26 M~ 3 ' 5 M can be 
as small as Mr < IOM0 for halo masses M h > 1.2 x 1O 7 M . 
However, thermal instability does not work in the haloes 
of smaller masses when Mr exceeds their baryon mass. As 
it seen from Fig. 7 this occurs at M h < 9 x 10 5 M Q . The 
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Figure 6. Jeans mass (upper plate) and temperature (lower plate) reached at one collision time in the postshocked gas as a function 
of the halo mass and the rcdshift when the haloes merge. White dotted line marks masses of the haloes whose compressed layers are 
gravitationally unstable. Scale bars in the right are for Jeans mass (upper) in solar units, and gas temperature (lower) in K. 



clouds formed under thermal instability are pressure sup- 
ported, but not necessarily gravitationally bound - it is seen 
from Fig. 6 that thermal instability forms gravitationally 
bound clouds of stellar masses only in mergings of massive 
haloes M £ (1 — 3) x 1O 7 M . Therefore not always thermal 
instability can give rise directly to star formation, however 
collisions of clouds formed through it can be a stimulating 
factor. 

The compressed layers gravitationally unstable against 
fragmentation, M > M^ r , lie in the halo mass range where 
HD cooling can keep gas isothermal with temperature close 
to Tomb- Therefore, once the fragments are formed, they 
can afterwards either collapse homogeneously or brake fur- 
ther on smaller masses as prescribed by the hierarchical 
fragmentation scenario (Hoyle 1953). The possibility of the 
fragmentation cascading down to lower scales still remains 
under debate (see, e.g., discussion in Coppi et al. 2001, 
Glover 2005). An essential condition for such a cascad- 
ing process is a pressure- free contraction regime. In prac- 
tice, gravitational collapse turns out to settle onto a slow 
quasi-hydrostatic regime, which results in an efficient sub- 



sonic damping of perturbations (Abel et al. 2002, Bromm 
ct al. 2002), and therefore in a strong suppress of fur- 
ther fragmentation. The mass of dense clumps which can 
be reached in the hierarchical fragmentation scenario gives, 
however, an estimate of the lower mass limit of protostellar 
condensations expected at given physical conditions - the 
opacity limit for fragmentation (Low & Lynden-Bell 1976, 
Rees 1976). This mass is determined by a maximum den- 
sity when the forming fragments become optically thick in 
the lines providing dominant radiative cooling. When HD 
molecules are the main cooling agent, a hierarchical frag- 
mentation becomes optically thick in HD rotational lines, 
i.e the optical depth of a fragment at line centre becomes 
thd = 1, at densities ~ 10 9 — 10 10 cm~ 3 . Even though at 
these densities most of hydrogen converts to molecular form 
in three-body collisions (Palla et al. 1983) its contribution 
to the total cooling remains negligible because of low tem- 
peratures. The corresponding Jeans mass at this stage is 
Mj ~ 3OT^j B n~ 1/2 M ~ 10" 3 (1 + zf /2 M Q (Vasiliev & 
Shchekinov 2005). For z — 10 — 20 this gives a relatively 
low mass limit: Mj ~ (0.03 — 0.1)M©. A similar estimate 
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Figure 7. Upper panel: The temperature dependence of the effective cooling function A(T) (dashed line), and its logarithmic derivative 
<E> = (d In A/9 In T) n (solid line). In a non-steady cooling medium temperature interval with <E> < 2 is thermally unstable. Initial density 
value is 1.5 cm~ 3 . Lower panel: Characteristic mass of condensations formed through thermal instability is shown by dotted lines. Solid 
lines show Jeans mass in the shocked layer. The solid and dotted lines on both panels are for halo masses M h = 9 x 10 5 ,10 6 ,2 x 10 6 ,3 x 
10 6 ,5 x 10 6 ,8 x 10 6 , 1.2 x 1O 7 M from top to bottom. Each curve is shown in the temperature interval from the initial value To to the 
lowest value reached at the end of compression t = t c for a given halo mass. Horizontal short-dashed lines show the baryon mass of the 
same set of halo masses in a reversed order - from bottom to top. 



for the mass of fragments formed behind the shock waves 
with velocities ^ 300 km s _1 from supernovae explosions 
was obtained by Uehara & Inutsuka (2000). One should 
stress though, that the condition thd = 1 implies that pho- 
tons from HD rotational lines experience in average only one 
scattering, which seems insufficient to lock radiation inside 
a contracting cloud and stop further decrease of Jeans mass. 
From this point of view, the condition thd =3 — 5 equivalent 
to about 10 scatterings of photons (Ivanov, 1973), gives ap- 
parently more accurate estimate of the opacity mass limit: 
the corresponding Jeans mass is Mj ~ (0.01 — 0.03)M Q . 
For comparison, the estimates of the opacity mass limit of 
the fragments determined by H2 cooling vary from 3Mq to 
about 1O 2 M0 (Nakamura & Umemura 1999) depending on 
the initial H2 abundance. Note, that the exact value of the 
opacity limit depends on dynamical regime of a collapsing 
cloud: the optical depth of a contracting cloud decreases 
with the velocity gradient as \dv/dr\~ (Sobolev, 1960). In 



numerical simulations \dv/dr\ ~ 2vt/Rj at late stages of the 
contraction (Abel et al. 2002), here Rj is the Jeans length. 
Abel et al. (2002) stop their similation when optical depth at 
line centre is 10. With accounting the factor 2 in the veloc- 
ity gradient this corresponds to the effective optical depth 
of 5 - nearly at the beginning of the stage when radiation 
becomes locked inside the contracting fragment. 



5 CONCLUSIONS 

In this paper we showed that 

(i) In the hierarchical scenario of galaxy formation merg- 
ings of massive dark matter haloes, M 8 x 10 6 [(1 + 
z)/20]~ 2 Mq develop dense baryon layers, where practically 
all the deuterium becomes confined into HD molecules, 
which then provide a very efficient radiative cooling; 
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(ii) Only slightly higher is the critical mass of the merg- 
ing haloes when the shocked baryon layers can be un- 
stable against gravitational fragmentation M > 10 7 [(1 + 
z)/20]~ 2 M e . 

(iii) In these conditions the shocked baryons can cool 
down to the lowest temperature T ~ Tcmb = 2.7(1 + z) 
K; 

(iv) The fragments formed through instability evolve then 
practically isothermally with T ~ Tcmb, and can in princi- 
ple undergo further sequencial fragmentation, which stops 
when the minimum mass Mq ~ 10 _3 (1 + z) 3//2 Mq is 
reached. 

(v) This minimum mass is smaller than the one expected 
when H2 molecules control radiative cooling (Nakamura & 
Umemura 2002, Ciardi & Ferrara 2005). This means that 
the first stars can be less massive in galaxies with masses 
M > 8x 10 6 [(l + z)/20]~ 2 M Q , where HD cooling dominates. 
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